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Abstract 

We review our recent studies on the Kondo effect in the tunneling phenomena 
through quantum dot systems. Numerical methods to calculate reliable tunneling 
conductance are developed. In the first place, a case in which electrons of odd 
number occupy the dot is studied, and experimental results are analyzed based on 
the calculated result. Tunneling anomaly in the even-number-electron occupation 
case, which is recently observed in experiment and is ascribed to the Kondo effect 
in the spin singlet-triplet cross over transition region, is also examined theoretically. 
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1 Introduction 



Quantum dot systems are now designed as the artificial magnetic impurity and 
are growing as a field of detailed experimental studies of the Kondo problems]!] 
. In this report we review our recent theoretical works on the Kondo effect in 
the tunneling phenomena through quantum dot [2-4]. 

After the papers pointed out the possibility to occur the Kondo effect in tun- 
neling through a quantum dot [5], many theoretical studies have been done 
on this problem[6]. The calculation of the tunneling conductance needs the 
dynamical excitation spectra. However we have not exact analytic calculation 
of the dynamical excitation spectra for the Kondo systems [7]. Two numerical 
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methods have been recently developed to calculate the tunneling conductance 
of the quantum dot systems. One is based on the numerical renormalization 
group technique (NRG) [8], and another is based on the Quantum Monte Carlo 
(QMC) technique [3]. Both techniques are known as reliable methods to cal- 
culate the dynamical excitation of the Kondo systems[7,9,10]. 

When the occupation number of electrons on the dot is odd, localized spin 
freedom appears on the dot, and it couples with the conduction electrons on 
the leads. At very low temperatures, we can expect the increase of the tun- 
neling conductance due to the resonance transmission via the Kondo peak in 
the density of states on the dot orbitals. This is the most typical example of 
the Kondo effect of the dot systems, and has been observed in many experi- 
mentsfl]. In the first part of this report, we present the theoretical calculation 
for this case[2] and compare it with experimental datafll]. Recently, anomaly 
in a region of an even electron number occupation case has been reported 
in experiment [12]. This phenomena is expected to relate to the Kondo effect 
in the spin crossover region of even occupation number case[4,13,14]. This 
problem is discussed in the second part of this report. 



2 Single Orbital Case 

We consider the following Hamiltonian, 



H = H i + H d + H i . d , (1) 

H i = J2 £ ^ C tka C vka, (2) 

(3) 

PO" <pcr,p'<j> 

Ht-d = -j= E {vpnkd+ a c^ ka + h.c.}. (4) 

V J- * puka 

The terms H e and H d represent the electron in the leads and the dot, respec- 
tively. The term gives the electron tunneling between the leads and the 
dot. The suffix ji = L(R) means the left(right) lead and dp means the dot 
orbital denoted by p. The quantity e^p corresponds to the energy of the or- 
bital, and it can be changed by applying gate voltage. The quantity U is the 
Coulomb interaction constant. 



At first we consider the most simplified model that the dot has a single or- 
bital. We abbreviate the suffix p. There will be many orbitals in dot in actual 
situations. Two orbitals case will be discussed in §3. We calculate the con- 
ductance, G, in the linear theory of the bias voltage. It is obtained from the 
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correlation function of the current operators [8]. But in the single orbital case, 
the calculation is reduced to the following expression [3], 

2e 2 f AT L T R df 

G = JTJ f7TT^ ( - ImGdd(e))( -^ )de ' (5) 

where Gdd(e) is Green's function of the dot orbital and /(e) is the Fermi 
distribution function. We use this formula in this section because the numerical 
calculation of G^e) is easier than the calculation of the current correlation 
function. Detailed comparison of both methods are made in Ref.[3]. Hereafter 
we assume that the leads have constant density of states from — D to D with 
D — 1. The hybridization strength is parameterized as T M = n\v^\ 2 p c , where 
p c = 1/2D is the density of states for the lead. 

In Fig.l, we show the conductance as a function of the gate voltage (e^) 
for various temperature cases [2]. This calculation is carried out by the NRG 
method. At higher temperatures we have paired Coulomb oscillation peaks 
at Ed = — U and 0. They grow without increase of their width, thus become 
very sharp as the temperature decreases. At the same time the peak positions 
shift slightly to Ed = —U/2 side. When the temperature decreases further, the 
intensity of the valley region between the two peaks gradually increases, and 
the peaks merge into a broad single peak at extremely low temperature. The 
characteristic temperature Tkm varies drastically as Ed changes. (We define 
Tkm following the usual definition, Tk = 4%(T = 0), for each Ed case). It is 
the lowest at the mid point of the two peaks, Ed = —U/2, and is denoted as 
T£ M , which is T KM /£7 = 5.34 x 10~ 7 for the parameter case used in Fig.l. 
The conductance at the mid point begin to increase at about 10T£ M , and the 
valley disappears at about 0.2T£ M . 

Goldhaber-Gordon et al. have shown the detailed temperature dependence of 
the conductance [11]. We directly compare the NRG data with experimental 
data in Fig. 2 [2]. The Coulomb repulsion has been estimated to be 1.9 meV. We 
have estimated T = 0.12meV =(1.3 x 10 3 mK), and thus T/nU = 1.9 x 1(T 2 
in the experimental situation. On horizontal axis the points Ed = 0, —U of 
the calculated conductance are set at V g = — 119, — 92meV, respectively, and 
factor 0.31 is multiplied to the NRG data to fit the experimental data at 
the lowest temperature. This factor 0.31 would be caused by the asymmetry 
Tl 7^ Tr. For the left hand side peak, the conductance data agrees very well 
with the experimental one in 100 mK ~ T ~ 1500mK. This agreement suggests 
that the behaviors in experiment are caused by the Kondo effect in the mid- 
temperature region shown in Fig.l. The Kondo temperature at the valley 
region seems to be less than 10 mK. 

There are several discrepant points. At high temperature region T > 2000mK, 
the conductance of the experimental data is larger than that of the NRG 
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Fig. 1. Conductance as a function of Ed for the parameter T/{ttU) = 1.0 x 10 

(u = i.o x io~ 2 ,r/^ = 1.0 x 10~ 4 ). 



EXP. NRG (r/U=0.0197t) 




V g (mV) 

Fig. 2. Comparison between the experimental and the calculated data. The ex- 
perimental data is reproduced from Fig. 2 of Ref.fll]. The parameter of the NRG 
calculation is chosen to be T/(ttU) = 1.9 x 10~ 2 
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Fig. 3. (a) Occupation numbers of the e and o orbitals at T = 0, (n e ), (n Q ), and 
the total occupation number on the dot, (b) Conductance at T = for the two 
channel case, Gf2- The symbols for 5 are common in (a) and (b). The dashed arrows 
show the increasing of 5. Parameters are Ju/U = —0.3, T c /U = T Q /U = 0.02. 

data. This might be caused by the multi-orbital effect. The conductance shows 
disagreement in the valley and the right hand peak positions. The change of 
the gate voltage on the dot will affect not only the potential of the dot, but 
also the hybridization strength between the dot and the lead states. 

The conductance at low temperature is strongly suppressed by application of 
the Zeeman field. The theoretical study has been carried out based on the 
QMC method [3], and the result agrees well with experimental one[f5]. 



3 Spin Crossover Transition 

Recently, Sasaki et al. observed low temperature tunneling anomaly in even- 
number-electron-occupation case[12]. They controlled the energy splitting of 
orbitals by tuning the magnetic field[16]. When the level splitting is gradually 
increased in the even electron number case, the energy of the spin triplet 
state increases compared with that of the singlet state. The electrons occupy 
different orbitals to gain Hund's coupling energy in the former state, while 
the electrons with opposite spin occupy the lower energy orbital in the latter 
state. It was suggested that the low temperature anomaly is related to the 
Kondo effect due to this spin crossover transition [13, 14]. However detailed 
calculation of the conductance has not been done. It is observed that a bump 
grows between the two Coulomb peaks, i.e., a third peak appears between the 
Coulomb peaks and it grows when the temperature decreases. This behavior is 
quite different from that of the odd number case discussed in previous section. 

We calculate the conductance around the singlet-triplet crossover region for 
a system with two orbitals[4]. The orbitals are denoted as even(p=e) and 
odd(p=o). The energy of e(o) orbital is defined as — 5/2 (sd + 5/2). The 
exchange term H cx = J H Ea^^^)^^)^^^^^ is added to 
the Hamiltonian Eq.(l). (J H < 0) 



5 




0.0 1 1 1 1 1 L 




(H-e<0/U 

Fig. 4. Conductance at various temperatures as a function of gate voltage in the 
two channel case. (a)6/U = 0.24, (b)6/U = 0.282, (c)S/U = 0.3. The symbols of 
the temperature in (b) and (c) are common. 

In Fig. 3 (a), we show the occupation numbers of the even and odd orbitals, 
(n c ) and (n ), on the dot. The parameters are Jr/U = —0.3, T c /U = T /U = 

0.027T. 

We assume the case in which both leads have two conduction channels, and 
each dot orbital connects to each conduction channel (two channel case). The 
conductance is given as the sum of ones from two channels. At very low tem- 
peratures it is given as, Gp2 = (2e 2 //i) X) P sm2 ( 7r (^p)/2), and is shown in 
Fig.3(b). 
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Fig. 5. Conductance as a function of gate voltage in the single channel case. Pa- 
rameters are the same to those of Fig.4(b). 



Here we define a quantity x = {^i — e^/U . In the ~ x ~ 1 region, an electron 
occupies the e-orbital, ((n c ),(n )) ~ (1,0). Next we see the 1 ~ x ~ 1.5 
region. At S/U = 0.24, the occupations on both orbitals are almost the same, 
((n e ), (n )) ~ (1, 1), due to the strong Hund's coupling. When S increases, the 
occupations in the 1 ~ x ~ 1.5 region gradually split to (2, 0) where the local 
spin state is a singlet. Therefore the conductance decreases from Ae 2 /h to 0. 
We stress that the splitting is suppressed around x — 1.5 as seen in Fig. 3 (a). 
For example at S/U = 0.28, the second electron tends to occupy the e-orbital 
when x sweeps to x — 1.2 from smaller x, then they redistribute to the e- and 
o-orbitals when x further sweeps to x — 1.5; (1,0) — > (1.6,0.3) — > (1.2,0.8). 
This redistribution reflects the stronger Hund's coupling near x = 1.5, at 
which the total occupation is close to two. The redistribution causes a bump 
in the conductance around x = 1.5 at low temperatures, as seen in 0.276 ~ 
S/U ~ 0.284 in Fig.3(b). 

In Fig. 4 we show the conductance at various temperatures around the lo- 
cal spin singlet-triplet degeneracy. The parameters are: (a) S/U = 0.24, (b), 
S/U = 0.282, (c) S/U = 0.3. There are the four Coulomb peaks on x = 
(fi — e$)/U ~ 0, 1, 2, 3 at high temperatures. When the temperature de- 
creases, the conductance in the region ~ x ~ 1 (2~x~3), where N d ~ 1 
(iV d ~ 3), increases to ~ 2e 2 /h, caused by the usual spin- 1/2 Kondo effect for 
the e-channel (o-channel). On the other hand, the region 1 ~ x ~ 2, where 
Nd ~ 2, is remarkable. At (a) S/U = 0.24, the conductance is large, about 
4e 2 /h. However, the large conductance appears at extremely low temperature 
(e.g., T/U = 10~ 8 corresponds to 10~ 3 mK for U ~ lOmeV). When S increases 
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to (b) S/U = 0.282, a "bump" emerges in the region 1.2 ~ x ~ 1.8 as the 
temperature decreases. For the case (c) S/U = 0.3, the conductance in the 
region 1 ~ x ~ 2 is nearly constant with a small value. The behaviors (a)-(c) 
are classified into, (a) the local spin triplet Kondo effect, (b) the singlet-triplet 
Kondo effect, and (c) the usual even-odd oscillations for the dot with large 
energy separation [1,8]. 

Next we consider a case in which leads have only one conduction channel, and 
the both dot orbitals connect to this single channel (single channel case). In 
such case one orbital will hybridize mainly with even combination of L and 
R lead states, and another mainly with odd combination of L and R states. 
The tunneling via different orbitals interferes in this case. At T = 0, the 
conductance is given as G-pi = (2e 2 //i) sin 2 (7r((n e ) — (n G ))/2) , and is small 
when (n e ) ~ (n ). Therefore the conductance will tend to zero as T decreases 
to extreme low temperature when Hund's rule coupling energy dominates. 
This is contrasted to the case (a) in Fig. 4. 

In the singlet-triplet Kondo effect case where the Kondo temperature is not 
low, the calculated temperature dependence of the conductance is not so dif- 
ferent from that of the two channel case(Fig.4(b)) as seen from Fig. 5. This 
is because (n e ) and (n ) have different and not-integer values as seen from 
Fig.3(a). 

As a summary, we have shown that the electron occupation on orbitals redis- 
tribute to gain Hund's coupling energy in the singlet-triplet crossover region 
when the potential deepens. This redistribution causes a bump, as seen in the 
experiment, in the conductance at low temperatures. 
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